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16  Abstract 

'The  object  of  the  research  was  to  investigate  the  amplification  of  a  Reynolds  number 
dependent  process  by  wave  distortion  and  the  possibility  of  applying  the  results  to  other 
similar  Revnolds  number  dependent  processes.  I  he  process  investigated  was  that  associated 
with  the  operation  of  a  constant  temperature  hot  wire  anemometer.  The  application  of  most 
interest  was  investigation  of  vaporization  limited  combustion,  the  type  combustion  typcially 
associated  with  liquid  propellant  rocket  engines. 


A  series  of  experiments  were  carried  out  to  determine  the  effect  of  wave  distortion  on 
a  Reynolds  number  dependent  process  and  to  establish  the  analogy  between  the  anemometer 
process  and  the  combustion  process.  Both  analytical  and  experimental  studies  were  involved, 
both  having  the  ultimate  objective  of  identifying  parametric  trends,  recognizing  behavior  common 
to  different  chamber  geometries  and  determining  stability  boundaries. 


The  results  indicate  a  high  degree  of  similarity  between  the  two  processes  and  the 
[X>ssibility  of  using  the  anemometer  system  to  investigate  combustion  instability.  The 
nonlinear  aspects  of  a  Reynolds  number  dependent  process  appear  to  be  the  dominant 
mechanisms  controlling  instability. _.'v 
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The  advice,  encouragement,  and  support  of  Richard  J.  Priem,  the  technical 
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Introduction  and  Background 


\ 

Thi-  object  of  this  report  is  to  summarize  the  research  and  the  signi¬ 
ficant  results  of  the  research  which  was  carried  out  undei  NASA  Grant 
NGR  43-003-015  from  January  1,  1972,  until  the  grant  termination  on 
Oc  t  obe t  31,  1979- 

The  purpose  of  the  research  was  to  study  high-frequency  combustion 
instability  in  liquid-fuel  rocket  eng i nes-- i  11s t ah i I i t y  in  which  the  com¬ 
bustion  process  is  coupled  to  acoustic  waves.  To  do  this,  new  experimenta 
and  analytical  tools  and  techniques  were  developed.  Basically,  it  was 
desired  to  develop  a  better  understanding  of  and  more  insight  into  the 
fundamental  phenomena  associated  with  combustion  instability,  with  the 
view  that  results  of  such  studies  would  not  be  confirmed  in  application 
to  liquid-fuel  combustor  instability,  (a  Reynolds  Number  dependent  pro¬ 
cess),  but  would  have  general  application  to  other  Reynolds  Number  depende 
processes.  These  other  processes  would  include  those  associated  with  all 
types  ot  combus ters-- so  1 i d-  as  well  as  liquid-fuel  rocket  engines,  jet 
engines,  and  furnaces. 

Combustion  instability  has  been  a  persistant  problem  in  both  liquid- 
and  solid- fuel  rocket-engine  development  programs.  In  liquid-fuel  rocket 
engines,  the  fuel  and  oxydizer  are  injected  into  the  combustion  chnmbei 
where  they  are  vaporized  and  combustion  takes  place.  The  resulting  hot 
gases  continue  through  the  chnmbei  and  are  accelerated  to  supersonic 
velocities  by  means  of  a  converging-diverging  nozzle  and  ultimately  leave 
the  engine.  This  process  is  not  smooth,  but  has  associated  with  it  back¬ 
ground  oscillations  which  are  detected  as  fluctuations  in  gasdvnamic 
variables  such  as  pressure,  density,  temperature,  and  so  forth.  (These 
oscillations  can  be  a  result  of  a  number  of  phenomena  such  as  turbulence 
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or  combustion  noise.)  The  combustion  process  and  the  oscillations  are 
coupled.  Any  oscillation  in  gasdynamic  variables  causes  an  oscillatory 
burning  rate;  and  the  oscillatory  burning  process  serves  as  the  energy 
source  to  sustain  gasdynamic  oscillations.  Undei  certain  conditions, 
when  the  energy  supplied  by  the  unsteady  combustion  process  is  sufficient 
to  overcome  the  energy  losses  in  t  he  gasdynam i c  field,  the  coupling  is 
se 1 f-sus ta i ned  and  the  oscillations  can  build  to  large  amplitudes,  a 
condition  referred  to  as  unstable  combustion.  So,  in  part,  the  question 
to  be  addressed  by  the  research  was,  why  and  how  the  random  oscillations 
which  normally  occur  during  operation  of  a  rocket  engine  change  character 
and  develop  into  an  organized  pattern  of  acoustic  oscillations  which  can 
build  in  magnitude  and  very  quickly  destroy  the  eng'ne  if  allowed  to 
con  t inue . 

When  this  research  was  initiated,  much  work  had  already  been  done 
experimental ly  and  analytically  in  the  area  of  combustion  instability, 
but  marked  difficulties  existed.  Direct  experimental  observation  of 
combustion-chamber  instabilities  was  very  costly  and  the  results  were 
difficult  to  interpret.  Experimental  modeling  techniques  were  not  entirely 
successful.  Some  things  were  known,  however.  For  example,  the  oscillations 
associated  with  unstable  combustion  usually  consisted  of  some  fundamental 
frequency  oscillation  distorted  by  an  oscillation  at  the  doubled,  or  second 
harmonic,  frequency.  Even  if  initially  the  instability  was  dominated  by 
a  single  frequency  oscillation,  it  would  quickly  become  harmonically  dis¬ 
torted  as  the  instability  developed.  Also,  the  oscillations  which  occurred 
could  be  most  closely  characterized,  in  frequency  and  form,  as  the  resonant 
oscillations  associated  with  a  closed-closed  cylindrical  chamber--the  funda¬ 
mental  oscillation  in  unstable  combustion  was  similar  to  the  first  tangential 
resonant  mode  of  a  closed-closed  cylindrical  chamber  of  appropriate  dimensions. 


In  addition,  some  initial  minimum  amplitude  of  gasdynamic  disturbance 

was  necessary  to  excite  instability. 

\ 

Analy t  ica  1  ly  ,  some  success'had  been  achieved  with  modeling.  It 
was  generally  accepted  that  combustion  instability  was  primarily  sensitive 
to  fluctuations  in  the  gas  pressure  and  the  velocity  of  the  fuel  drops 
relative  to  the  gas.  The  problem  of  pressure-sensitive  combustion  in¬ 
stability  had  a  long  history  of  successful  analytical  treatment;  but, 
because  of  its  greater  mathematical  complexity,  much  less  work  had  been 
done  with  analysis  of  velocity-sensitive  combustion  instability. 

Acquiring  the  understanding  of  combustion  instability  needed  to 
formulate  rational  design  procedures  depended  critically  on  the  develop 
merit  of  realistic  modeling  techniques,  and  it  was  toward  this  end  that 
this  research  was  begun.  Initially  the  emphasis  was  on  development  ot  an 
experimental  analog  of  the  combustion-acoustics  interaction.  Later,  ana¬ 
lytical  techniques  were  developed,  partially  in  suppoi  t  ot  the  exper  imental 
work  and  partially  as  a  valuable  study  in  its  own  right. 

Several  specific  observations  led  to  the  initiation  of  this  research. 
First  Priem  [l]  established  a  model  in  which  combustion  is  viewed  as 
vaporizat ion-1 imi ted  (i.e.,  as  a  velocity-sensitive  or  Reynolds  number 
dependent  process^.  Later,  He i dman  [2,  3]  showed  that  distorted  acoustic 
oscillations  (fundamental  frequency  oscillations  plus  high  harmonic 
"distortion")  affect  the  open-loop  response  of  this  vapor i zat i on- 1 imi ted 
combustion.  The  response  is  a  measure  of  the  degree  of  reinforcement  of 
the  acoustic  oscillations  and  the  combustion  process.  It  is  important  to 
know  how  the  i ns tantaneous  gasdynamic  acoustic  field  combines  with  the 
acoustic  field  resulting  from  the  combustion  process.  One  of  the  better 
ways  of  relating  these  two  fields  is  to  find  the  correlation  between  the 
pressures  of  each.  Such  a  correlation  could  then  be  used  to  determine  the 
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closed-loop  response.  But  no  explicit  expression  for  the  acoustic  pres¬ 
sure  which  results  from  a  particular  combustion  process  was  available; 
hence.  Heidmann  analytically  compared  the  existing  pressure  with  the 
burning  rate  of  the  combustion  process,  i.e.,  he  formed  an  open- loop 
response.  He  used  harmonic  response  factors  to  indicate  the  degree  of 
reinforcement  of  the  acoustic  pressure  by  the  combustion  process  and 
found  that  the  response  factors  depended  upon  the  amplitudes  of  the  har¬ 
monic  components  and  the  phase  between  the  acoustic  pressure  and  burning 
rate  harmonic  components.  Lastly,  Hribar  [A]  noticed  the  similarity  of 
the  open-loop  responses  of  a  liquid-fuel  droplet's  vaporization  rate  and 
a  cons tant- temperature  hot-wire  anemometer's  energy  transfer  rate  to  the 
same  gasdynamic  environment. 

Combining  the  above  observations,  Hribar  [I*]  proposed  the  development 
of  an  analog  technique  for  studying  combustion  instability,  in  which  a 
cons  tan t - tempo  rat ure  hot-wire  anemometer  would  servo  as  the  analog  to 
combustion.  The  advantages  to  such  a  system  would  be  many.  The  analog 
would  not  involve  an  actual  combustion  process  so  that  the  many  problems 
associated  with  making  measurements  in  a  hot,  burning  mixture  would  be 
avoided.  The  system  could  be  more  easily  controlled;  studies  could  be 
done  with  the  analog  which  would  not  be  possible  using  an  actual  combustion 
situation.  (For  example,  parameters  could  be  varied  independently,  a 
feature  not  possible  in  an  actual  combustion  situation.)  Also,  the  analog 
technique  would  be  less  expensive  than  actual  experimental  combustion  studie 

Because  of  the  above,  the  initial  objective  of  the  research  was  to 
develop  an  analog  system  for  the  study  of  combustion  instabilities  in  liquid 
propellant  rocket  engines,  and  then  to  use  the  analog  to  investigate  the 
basic  mechanisms  governing  these  instabilities.  Later,  the  research  was 
expanded  to  include  mathematical  analysis  techniques,  initially  to  determine 


whether  or  not  certain  effects  that  were  observed  experimentally  with  the 

analog  could  be  predicted  mathemat ical ly.  These  analytical  studies  re- 

v 

suited  in  improved  insight  i nto'ana ly t i ca 1  modeling  techniques. 

What  follows  in  this  report  is  a  chronol og i ca I  description  of  the 
research  and  the  significant  results.  Most  of  what  was  done  has  been 
documented  in  various  papers,  theses,  and  dissertations.  Hence,  specific 
details  of  the  research  will  be  kept  to  a  minimum  and  the  appropriate 
references  will  be  mentioned  in  which  more  detailed  descriptions  can  be 
found . 
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Research 


As  was  stated  earlier,  the  initial  objective  of  the  research  was 
to  develop  an  analog  system  for  the  study  of  high  frequency  combustion 
instabilities  in  liquid-propellant  rocket  combustion  chambers.  Be.fore 
considering  the  analog,  an  understanding  of  the  phenomena  being  modeled 


is  necessary. 

Combustion  Mode  I .  To  characterize  the  combustion  process,  the 
combustion  model  due  to  Priem  [l],  in  which  the  combustion  is  viewed  as 
vapor i za t i on- 1 imi ted ,  was  used.  In  this  model  the  burning  rate  is  equal 
to  the  vaporization  rate  of  the  liquid  propellant  droplets,  and  this  vapo¬ 
rization  or  burning  rate  W  is  given  as 

W  =  Ci  +  C ■  Re  '  ’  (la) 


where 


and 


2  PM  V  ,  p 

R~f  p  ,n  p^p- 

U  V 

0.6  Sc'  *  ''  :tDM,PC  4  p 

R  T  P  n  P  -  P 

u  v 


pD  |V  -  V  ,| 
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where  D  is  the  inside  diameter  of  the  combustion  chamber,  M.  is  the  mole- 

cular  weight  of  the  propellant,  V  is  the  molecular  diffusion  coefficient, 

C,  is  the  concentrat ion  of  propellant  drops,  R  is  the  universal  gas  con- 
A  u 

stant,  T  is  the  gas  temperature,  P  is  the  gas  pressure,  P^  is  the  vapor 
pressure  of  the  fuel,  S  is  the  Schmidt  number,  p  is  the  gas  density, 

V  is  the  gas  velocity,  V,  is  the  propellant  droplet  velocity  and  u  is 
the  dynamic  viscosity.  Equations  (1)  indicate  that  the  vaporization  rate 
is  convection  limited  (Reynolds  number  dependent)  with  additional  sensitivity 
to  the  gasdynamic  and  droplet  variables. 
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In  Priem's  model,  the  burning  rate  is  sensitive  to  the  magnitude  of 
the  surrounding  gas  velocity  with  respect  to  the  drop  velocity,  i.e.,  it 
is  sensitive  to  a  "rectified"  relative  velocity.  This  causes  higher 
harmonic  frequency  components  of  the  burning  rate  to  occur  in  response 
to  fundamental  frequency  perturbations  of  the  gasdynamic  field.  In  addition 
to  this  method  of  generating  higher  harmonic  components  of  the  burning  rate, 
and  consequently  of  the  gasdynamic  variables,  nonlinear  gasdynamic  effects 
will  also  generate  higher  harmonic  components  of  the  gasdynamic  variables. 

As  was  mentioned  earlier,  during  combustion  in  liquid-propellant  rocket 
engines,  the  combustion  process  and  the  gasdynamic  processes  are  coupled. 

An\  oscillation  in  gasdynamic  variables  causes  an  oscillatory  burninq  rate; 
and  the  oscillatory  burning  process  serves  as  the  energy  source  to  sustain 
gas  dynamic  oscillations.  Under  certain  conditions,  when  the  energy  supplied 
bv  an  unsteady  combustion  process  is  sufficient  to  overcome  the  energy  losses 
in  the  gasdynamic  field,  the  coupling  is  se I f- sus ta i ned  and  the  oscillations 
can  build  to  large  amplitudes,  a  condition  referred  to  as  unstable  combustion. 

The  Analog.  The  analog  selected  to  model  the  vapor i zat i on- 1 i mi  ted 
combustion  process  was  a  hot-wire  sensor  and  the  associated  series  output 
of  a  constant  temperature  hot-wire  anemometer,  a  phase  change  or  time  delay 
device,  a  power  amplifier,  and  an  acoustic  driver.  The  heart  of  the  analog 
is  the  anemometer  which  applies  a  voltage  across  the  hot-wire,  causing  it 
to  heat  to  some  preset  temperature  (resistance)  level.  Anv  flow  which  passes 
over  the  wire  tends  to  cool  it,  but  the  anemometer  responds  to  counteract 
this  by  increasing  the  voltage  across  the  wire.  The  bridge  voltage  of  the 
anemometer  E  is  related  to  the  flow  over  the  wire.  This  relationship  can  be 
expressed  as 

c  “  - ^ -  LV.3  -r  c4  Ke^ 

w 


+  RS>' 
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who  r  o 


nL  kf  [T,/T] 

w  f  t 

-  T]A 

(2b) 

iL  k  [T  ,T] 

w  f  t  / 

-1'  fTw 

-  T]B 

(2c) 

and 

PfD  V 

Re  (  =  — -  (2d) 

f  uf 

where  D  is  the  diameter  of  the  hot-wire,  R  is  the  operatinq  resistance 
w  w 

of  the  hot-wire, R  ..  is  the  resistance  of  the  cable  connecting  the  hot-wire 
cb  I 

and  probe  to  the  anemometer,  R^  is  the  probe  resistance,  R  is  the  appropriate 
anemone te r  resistance,  m.  A,  and  B,  are  coefficients  whose  exact  values  are 
determined  by  calibration  of  the  individual  hot-wire  (typical  values  are 
m  =  A  =  0.5205  and  B  =  0.6729;  details  on  hot-wire  calibration  can 

be  found  in  [5],  [6],  [7],  [8],  L  is  the  length  of  the  hot-wire,  is 

W  I 

the  gas  (air)  thermal  conductivity  evaluated  at  the  film  temperature  T  = 

[T  +  T^]/2,  T  is  the  gas  temperature,  T  is  the  hot-wire  operating  tempera¬ 
ture,  .  .  is  the  gas  density  evaluated  at  the  film  temperature,  V  is  the 
component  of  the  gas  velocity  normal  to  the  stationary  hot-wire  and  g ,  is 
the  dynamic  viscosity  evaluated  at  the  film  temperature.  (Discussions  of 
the  development  of  Equations  (2)  are  given  in  several  references  by  Purdv, 
Ventrice,  and  Fang  [5,  6,  8,  9].)  Equations  (2)  indicate  that  the  square 
of  the  anemometer  output  voltage  is  analogous  to  the  convection  limited 
burning  rate  W  given  by  Equations  (I).  To  complete  the  analogy,  a  phase 
change  or  time  delay  device  is  used  as  the  analog  of  the  phase  change  or 
tine  delay  associated  with  combustion,  a  power  amplifier  is  used  as  the 
analogof  the  energy  released  during  the  burning  of  the  propellant,  and 
an  acoustic  driver  is  used  to  convert  the  energy  associated  with  the 
anemometer  and  power  amplifier  to  gasdynamic  energy. 
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In  the  analog  model,  the  anemometer  output  is  dependent  upon  the 

magnitude  of  the  gas  velocity  normal  to  the  hot-wire,  i.e.,  it  is  sensi- 

\ 

tive  to  a  "rectified"  relative  velocity.  As  in  combustion,  this  causes 
higher  harmonic  frequency  components  of  the  anemometer  output  to  occur 
in  response  to  fundamental  frequency  oscillations  in  the  gas  velocity 
field.  In  addition  to  the  nonlinearity  of  the  analog  process  resulting 
in  higher  harmonic  components  of  the  anemometer  output  voltage,  and  con¬ 
sequently  of  the  gasdynamic  variables,  nonlinear  gasdynamic  effects  also 
generate  higher  harmonic  perturbat i ons  of  the  gasdynamic  variables. 

A  series  of  events,  analogous  to  those  which  occur  in  the  combustion 

system,  occur  in  the  analog  system.  Acoustic  perturbations  in  the  gas¬ 

dynamic  field  are  sensed  by  the  hot-wire.  The  resulting  anemometer  output 
signal  is  amplified,  after  phase  change  or  time  delay,  and  fed  to  the 
acoustic  driver,  which  generates  gasdynamic  perturbations  proportional 
to  the  signal  to  it.  These  new  acoustic  perturbations  combine  with  the 

previous  ones  to  form  a  new  gasdynamic  field.  As  in  combustion,  in  some 

situations  the  perturbations  will  reinforce  each  other  and  be  sustained 
or  increased  in  intensity;  in  other  cases  they  will  die  out. 

A  comment  needs  to  be  made  concerning  the  anemometer  output  signal 
mentioned  above.  The  anemometer  actually  has  two  output  signals--E,  the 
bridge  voltage  and  E ' ,  the  time-dependent  or  fluctuating  component  of  the 
bridge  voltage,  where  E  =  E  +  E'.  (The  component  E  is  the  time-mean 
component.)  It  is  E2  which  is  analogous  to  the  burning  rate  W  and  E2 
also  has  a  time-mean  and  a  fluctuating  component,  i.e.,  E2  =  (1^)  +  (E2) 1  . 

(The  bar  over  a  property  is  used  to  indicate  the  time-mean  component  and 
the  prime  designates  the  time-dependent  component.) 

From  the  above  it  can  be  seen  that  E2  =  E2  +  2EE'  +  (E1)2  =  (E2)  +  (E2)  '  . 
Only  the  fluctuating  component  of  E 2  would  drive  acoustic  oscillations; 


9 


hence,  ( E ‘ )  '  would  be  the  required  siqnal  to  be  used  in  the  analoq  system. 
However,  when  E'  is  small  relative  to  F  (which  is  the  case  for  the  analog 
system).  Fang  [8]  has  shown  that  E  1  character  i  zes  (£•  )  1  to  a  satisfactory 
degree  with  respect  to  amplitudes,  frequencies,  and  phases.  Hence,  E1 
was  used  as  the  feedback  signal  in  the  loop. 

Initial  Research  Plan.  At  the  outset  of  the  research,  the  above  des¬ 
cription  of  the  analog  had  not  been  fully  conceptualized  and  the  validity 
of  the  description  had  not  been  established.  As  the  research  progressed, 
the  above  description  evolved. 

To  initiate  the  analog  studies,  the  following  was  proposed. 

1.  To  analytically  perform  an  open-loop  response  factor  analysis 
of  the  analog  to  see  if  the  same  results  would  be  obtained  as 
Heidmann  obtained  for  vapor i za t i on- I i mi  ted  combustion.  For 
these  studies  the  pressure  associated  with  an  acoustic  field 
would  be  compared  with  the  resulting  anemometer  output  (analogous 
to  the  burning  rate  in  combustion)  and  a  response  factor  ob¬ 
tained.  Various  amounts  and  types  of  distortion  would  be 

added  to  determine  stability  predictions.  This  process  would 
establish  similarity  between  combustion  and  the  analog. 

2.  If  the  results  of  the  above  were  in  agreement  with  Heidmann, 
to  experimentally  perform  an  open-loop  analysis  of  the  ana¬ 
log  to  determine  if  the  analog  would  perform  as  predicted 

by  theory. 


3.  To  operate  the  analog  system  in  the  closed-loop  mode  to 
determine  if  the  open-loop  predictions  of  closed-loop  be¬ 
havior  were  valid.  The  above  three  items,  if  successful, 
would  confirm  the  validity  of  Heidmann's  analytical  technique 
and  the  prediction  that  distorted  acoustic  waves  would  drive 
combustion,  or  the  analog  process,  to  instability.  Once  all 
of  the  above  was  accomplished,  the  analog  system  could  then 
be  used  for  further  experiments  designed  to  identify  other 
factors  important  to  unstable  combustion. 


r t i cal  Open-Loop  Studies.  The  nonlinear,  in-phase  response  factor 


developed  and  used  by  Heidmann  is  defined  as 

r"  W*  p1  d  (ut) 

/  ‘  1  (p1  )  2  d  (ait) 
o 

where  W  is  the  time-dependent,  dimensionless  perturbation  in  the  normalized 


burning  rate  (W  W'/W)  and  p1  is  the  time-dependent,  dimensionless 
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comb i not i ons  of  harmonic  amplitude  and  phase  angle  would  result  in 

instability.  The  pure  sinusoidal  case  was  always  predicted  to  be  stable. 

v 

More  detailed  discussions  of  the  analytical  open- loop  studies  can  be 
found  in  [b],  [6 "| .  [9]  and  [  11  ] . 

Experimental  Open-loop  Studies.  Since  the  analytical  open-loop 
studies  were  successful,  the  experimental  open-loop  studies  were  carried 
out.  The  experimental  technique  for  doing  this  was  involved  and  is  dis¬ 
cussed  in  detail  in  [S]  and  [6], 

A  simplified  schematic  diagram  of  the  experimental  system  used  is 
shown  in  Figure  1.  The  object  was  to  correlate  an  electrical  signal, 
which  was  proportional  to  the  fluctuations  in  the  pressure  of  the  acoustic 
field  in  the  vicinity  of  the  hot-wire,  with  the  fluctuations  in  the  anemo¬ 
meter  output  voltage  and,  from  this,  to  calculate  the  response  factor  for 
that  particular  acoustic  field.  A  variety  of  acoustic  fields  were  used 
starting  with  an  undistorted  field  at  frequency  f^  and  gradually  adding 
to  this  distortion  at  2f  _.  so  as  to  parallel  the  work  done  ana  1 y t i ca I  I y . 

To  obtain  response  factors  the  system  illustrated  in  the  data  acqui¬ 
sition  system  of  Figure  I  was  used.  No  instrument  was  available  to  directly 
determine  a  response  factor  so  a  response  factor  was  calculated  from  the 
correlation  coefficient  R.  Again,  details  of  this  can  be  found  in  ['>]  and 
[6]. 

The  results  of  this  aspect  of  the  research  were  important  in  estab¬ 
lishing  the  analog  as  a  potential  research  tool .  First  of  all,  the  experi¬ 
ments  showed  that  it  was  possible  to  set  up  a  physical  system  which  could 
generate  the  same  acoustic  environment  for  a  hot-wire  probe  that  was  in¬ 
vestigated  analytically  and  which  could  be  used  to  determine  the  response 
of  the  hot-wire  voltage  to  that  environment.  Secondly,  the  experimental 
results  were  qua  1 i tat i ve ly  the  same  as  the  analytical  resu 1 1 s-- the  response 
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of  the  system  to  pure  sinusoidal  acoustic  vibration  was  small,  even  when 
the  magnitude  of  the  acoustic  pressure  was  large;  the  response  could  be 
increased  by  as  much  as  an  order  of*  magnitude  with  respect  to  the  sinu¬ 
soidal  case  by  the  addition  of  distortion;  and  the  amplitude  and  phase 
of  the  distortion  component,  relative  to  the  fundamental  component,  were 
the  dominant  factors  in  the  increase  in  the  response. 

Some  problems  did  exist.  The  distorted  acoustic  fields  could  not 
be  totally  character i zed  by  the  equations  resulting  from  the  solution  of 
the  inviscid  linear  wave  equation.  Viscous  effects  should  probably  have 
been  included  in  the  analysis  to  account  for  acoustic  streaming  and  any 
other  secondary  phenomena  which  might  have  been  occurring.  Also,  there 
might  have  been  some  problems  associated  with  the  assumption  used  to  obtain 
calibration  constants  for  the  hot-wire.  The  assumption  was  that  the  steady- 
state,  steady-flow  calibration  of  a  hot-wire  is  applicable  even  when  the 
wire  is  used  in  an  unsteady  situation. 

Other  discussions  of  the  experimental  open-loop  investigations  can  be 
found  in  [5],  [6],  and  [ll]. 

Experimental  Closed-Loop  Studies.  The  closed- loop  experimental  studies 
posed  many  difficulties.  Figure  2  is  a  very  simplified  schematic  diagram 
of  the  experimental  system.  The  test  section's  length-to-diameter  ratio 
could  be  adjusted  to  two  different  values--one  in  which  twice  the  frequency 
of  the  first  tangential  resonant  mode  of  vibration  was  equal  to  the  frequency 
of  the  second  tangent ia l-second  longitudinal  resonant  mode  of  vibration, 

2f]T  =  f,_  ;  and  one  in  which  twice  the  frequency  ot  the  first  tangential 

was  not  equal  to  a  resonant  frequency  of  the  chamber,  2f  _  =  f  .  Two  chamber 
were  used  to  establish  the  effect  of  chamber  dimensions  on  instability. 

Referring  to  Figure  2,  the  gas-filled  chamber  could  be  excited  by  moving 
the  switch  into  the  upper  position  so  that  the  function  generator/ampl i f ier 


Figure  2.  Simplified  Schematic  Diagram  of  the  Closed-Loop 
E  xpe  r  i  men  t  a  1  System 


system  supplied  the  excitation  signal  tv'  the  driver  ^open-loop  operat  i  on)  . 
By  switching  tv'  its  I  owe  r  position,  c  lv'sv'vi- loop  operation  could  bo  obtained 

One  objective  ot  the  closed- loop  invest i gat  ions  was  tv'  confirm  the 
results  of  the  open- loop  work,  i.e,.  if  the  system  was  operating  with  a 
certain  pattern  of  distorted  acoustic  osc i 1 lat Ions ,  would  it  hecv'me  unstably 
This  was  not  the  only  objective,  however.  Anothoi  objective  was  tv'  allvxv 
the  system  tv’  function  over  a  wide  range  of  variables  and  determine  under 
what  conditions  instability  developed  and  tv'  determine  it  these  conditions 
were  consistent  with  those  associated  with  unstable  combustion  and ,  from 
this,  tv'  perhaps  gain  insight  intv'  why  these  particular  conditions  develop. 

Tv'  obtain  a  better  understanding  ot  closed- loop  operation,  the  system 
was  initially  operated  using  a  linear  feedback  device--the  output  ot  a 
microphone.  This  was  a  pressure-sensitive  linear  feedback  process  instead 
of  the  velocity-sensitive  nonlinear  feedback  process  of  the  anemometer 
analog  system.  Using  the  microphone,  the  charac t er i s t i cs  ot  vlosed  loop 
system  behavior  at  a  single  frequency  could  he  studied  and,  ultimately,  a 
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comparison  could  ho  mode  of  pressure  and  velocity  sensitive  feedback 

processes  and  the  pressure  and  velocity  sensitivity  of  the  combustion 

\ 

process.  \ 

A.  Pressure-Sens i t i ve  C losed- Loop  L xpe r i men t  s 

Referring  to  Figure  2,  the  general  procedure  followed  was  to  excite 
the  chamber  in  the  open-loop  mode  with  some  predetermined  signal  by  means 
of  the  function  generator /amp  1 i f i e r  system.  The  switch  was  then  moved 
to  its  alternate  (closed- loop)  operating  position  and  the  system  was 
observed  to  see  what  would  evolve.  Three  basic  responses  could  develop: 
the  acoustic  pressure  in  the  chamber  could  die  out;  the  acoustic  pressure 
could  build  in  intensity  to  such  a  level  that  the  switch  had  to  be  moved 
back  to  its  open- loop  position  to  avoid  burning  out  the  fuse  in  the  driver 
circuit;  or,  finally,  some  type  of  steady-state  acoustic  pressure  field 
could  be  attained.  It  was  the  conditions  under  which  the  third  response 
occurred  that  were  of  interest.  This  represented  a  neutral  stability 
situation.  During  the  course  ot  the  experiments,  several  parameters  were 
altered  sys temat i ca 1 1 y  to  determine  their  effect  on  the  closed-loop  operation 
ot  the  system.  These  were  the  open-loop  signal  from  the  function  generator/ 
amplifier  system  to  the  driver,  the  amount  of  time  delay  or  amount  of  phase 
shift,  the  amount  of  amp  1 i f i cat i on  or  gain,  and  the  length-to-diameter  ratio 
of  the  chamber. 

The  overall  results  of  the  linear  feedback  study  were  as  follows: 

1.  The  c 1 osed- I oop  sound  field  was  independent  of  the  character  of 

the  initial  exciting  sound  field.  The  same  closed-loop  sound  field  resulted 
for  a  particular  amount  of  gain  and  time  delay  or  phase  change  in  the  loop, 
no  matter  what  the  chat acterist ics  of  the  initial  open-loop  sound  field. 

2.  For  a  particular  sound  field  to  be  sel f-sus tained  in  closed-loop 
operation,  the  least  amp  I i f i ca t i on  was  required  when  the  existing  and 
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feedback  sound  fields  were  in  phase  with  each  other.  If  these  two  fields 
were  not  in  phase  with  each  other,  more  than  this  minimum  feedback- loop 
amp  1 i f i c a t i on  was  required  in  order  to  have  se 1 f-sus tenance-- the  greater 
the  phase  difference  between  them,  the  greater  the  required  feedback  gain. 

3.  In  general,  the  higher  the  gain,  the  higher  the  level  of  the  closed- 
loop  soumi  field.  If  se 1 f-sus tenance  did  not  occur  at  a  certain  phase  dif¬ 
ference  between  the  existing  and  feedback  field,  increasing  the  gain  or 
adjusting  the  phase  difference  so  that  the  overall  phase  change  through  the 
loop  was  closer  to  a  net  zero  would  make  the  system  self-sustaining. 

4.  Coexistence  of  two  or  more  frequencies  was  possible  whenever  the 
gain  was  sufficient  for  each  of  them  to  be  se 1 f-sus ta i ned.  They  were  not 
normally  integral  multiples  of  each  other. 

5.  None  but  resonant  frequencies  were  self-sustaining. 

6.  The  length- to-diameter  ratio  of  the  chamber  did  not  effect  the 
character  of  the  results. 

P .  Ve  1  oc i t v-Sens i t i ve  C 1 osed- Loop  E xpe r i me n t s 

The  general  procedure  followed  was  similar  to  that  of  the  linear 
expe r i men t s - - the  chamber  was  excited  in  t he  open-loop  mode  with  some  pre¬ 
determined  signal  by  means  of  the  function  gene ra tor/amp I i f i er  system  (see 
Figure  2).  The  switch  was  then  moved  to  its  closed-loop  position  and  the 
system  was  observed  to  determine  what  type  sound  field  would  evolve. 

The  same  three  kinds  of  response  could  develop:  the  sound  field 
could  die  out;  the  sound  field  could  continuously  build  in  intensity;  or 
it  could  reach  some  steady-state  operating  level.  The  same  parameters 
could  be  altered--the  open-loop  excitation  signal,  the  time  delay  or  phase 
shift,  the  gain,  and  the  length-to-diameter  ratio  of  the  chamber. 

It  needs  to  be  recalled  here  that  the  anemometer 1 s  response  is  non¬ 
linear— it  is  related  to  a  "rectified"  relative  velocity.  Because  of  this, 


its  output  tends  to  be  dominated  by  a  component  at  2f  when  exposed  to  a 
velocity  field  of  frequency  f,  i.e.,  if  the  chamber  was  initially  excited 
with  an  '.T  sound  field  of  frequency  f ,  T  ,  the  anemometer  would  have  as  its 
output  a  siqnal  with  a  stronq  2f  ^  component. 

It  was  found  that  three  types  of  so  I  f -sus ta i ned  operation  could  occur 
in  the  chamber  for  which  2f  =  f,^.  depending  on  the  time  delay  or 

phase  change  in  the  loop:  (l)  steady  in  waveform  and  level,  having  fre¬ 
quency  components  at  f  and  at  its  higher  harmonics;  (2)  unsteady  in 

waveform  and  level,  having  frequency  components  at  f  and  its  higher 

1 T 

harmonics;  and  (3)  transitional  in  waveform  and  level  having  frequency 
components  at  f , ^  f , ^  and  higher  harmonics  of  each.  Figure  3  is  a 
map  of  the  1  owe  r  limit  of  gain  for  self-sustenance  at  various  values  of 
time  delay,  with  the  type  of  sound  field  resulting  indicated  by  the 
symbols.  These  results  can  be  partially  explained  using  the  results 
of  the  1 i near- feedback  experiments  and  calculations  of  the  net  phase 
change  through  the  loop  for  the  fundamental  and  second  harmonic  frequencies, 
the  results  of  which  are  listed  in  Table  I. 


Table  I.  Net  Phase  Change  Through  the  Loop  at 
Various  Time  Delays  for  f^  and  2f 
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Figure  3-  Se  1  f - Sus ta i ned  Map  of  the  Hot-Wire  Closed-Loop 

System  ( 2 f  T  =  f  T  ,  chamber) 

IT  .’T,  \l 

First,  Figure  3  follows  a  repetitive  pattern.  The  basic  pattern  is 
that  recorded  between  22.69  and  23.70  msec.  This  span  of  time  delay, 

1.01  msec,  is  equal  to  one  period  of  a  signal  at  f  ^..  Starting  at  22.69 
msec,  the  net  phase  changes  for  the  f  and  2f^  components  are  both  359 
degrees.  Since  the  existing  and  feedback  sound  fields  are  almost  in  phase, 
a  minimum  amount  of  gain  is  needed  in  the  loop.  At  22.80  msec,  where  the 
net  phase  changes  are  39  anil  78  degrees  respectively,  more  gain  is  needed 
for  self-sustenance.  At  23-21  msec,  the  net  phase  changes  are  1 83  and  7, 
respectively.  Since  the  fundamental  feedback  component  is  almost  280 
degrees  out  of  phase  with  the  existing  fundamental  component,  the  gain 
required  becomes  a  maximum.  When  the  existing  and  feedback  sound  fields  are 
180  degrees  out  of  phase,  the  level  of  the  feedback  sound  field  must  be 
greater  than  that  of  the  existing  field  if  self-sustenance  is  to  be  attained 
the  sound  fields  at  f  are  constantly  interfering  with  each  other.  The 
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level  of  the  se 1 f -sus ta i ned  sound  field  is  unsteady  but  the  frequency 

components  remain  the  same.  This  behavior  resulted  in  the  closed-loop 

v 

sound  field  appearing  "unsteady."  At  a  time  delay  ot  23.3^  msec--another 
relative  minimum  reqion--the  net  phase  change  for  the  f , j  and  2f,  com¬ 
ponents  are  2  3 2  and  !0;t  degrees,  respectively.  Both  values  of  time  delay 
are  s i qn i f i c an t  1  y  different  than  1  SO  degrees  so  less  gain  is  needed  for 
self-sustenance,  but  not  as  little  gain  as  at  23-70  msec  where  the  phase 
chan ges  have  returned  tv'  360  and  359  degrees,  respectively. 

One  time  delay  in  Table  1  has  not  been  men t i oned--23 . 00  msec,  at  which 
the  net  phase  changes  are  111  and  222  degrees,  respectively.  These  values 
of  net  phase  change  differ  from  zero  bv  more  than  they  do  at  22.80  msec 
and,  hence,  should  require  more  gain  instead  of  less  for  se 1 f -sustenance . 
However ,  another  factor  effects  the  results  at  22.80  msec.  As  mentioned 
earlier,  a  '.T,  ;  L.  sound  field  and  its  higher  harmonics  sound  fields  existed 
along  with  the  ;T  sound  field  and  its  higher  harmonic  sound  fields  in  this 
transitional  region.  This  second  series  of  sound  fields  reduces  the  gain 
needed  for  self-sustenance.  This  new  sound  field  appeared  in  some  cases 
due  to  the  following  character  i  s  t  ics :  (l)  f,^  ^  was  1135  Hz  while  f 

was  2262  Hz.  Hence,  f  was  approximate ly  twice  f,_  .  This  is  similar 

.'I  1  I  ,  » L 

to  2f^  =  f  -j-  .  It  was  concluded  that  whenever  a  resonant  frequency  was 

twice  another  resonant  frequency,  a  self-sustained  closed-loop  sound  field 
would  be  possible  if  the  gain  was  high  enough  and  the  time  delay  was  close 
enough  to  an  optimal  value.  For  this  reason,  the  appearance  ot  t  ^  ^  and 
its  higher  harmonics,  at  certain  time  delays,  was  to  be  expected.  (2)  The 
IT,  1L  sound  field  requires  more  energy  to  drive  (it  is  less  responsive) 
than  a  IT  sound  field.  Also,  f  is  not  exactly  twice  f;T  Hence, 

for  most  time  delays,  the  self-sustained  sound  field  is  naturally  dominated 


bv  f 


and  its  harmonic  harmonics  and  t 


does  not  develop. 


Another  significant  finding  of  the  nonlinear  experiments  was  that 
the  resultant  closed-loop  sound  fields  do  not  depend  on  the  form  (fre¬ 
quencies  and  phases)  of  the  exciting  sound  field.  However,  the  gain 
required  to  initiate  self-sustenance  did  vary  with  the  waveform  and 
level  of  the  exciting  sound  field.  When  the  open-loop  sound-pressure 
levels  for  various  resonant  modes  were  the  same,  the  minimum  gain  required 
to  initiate  se 1 f- sus tenance  occurred  when  the  exciting  sound  field  was 
the  T  field.  Also,  for  a  given  exciting  sound  field,  the  higher  the 
sound-pressure  level,  the  lower  the  gain  required  to  initiate  self-sustenance. 

Other  closed-loop  nonlinear  feedback  experiments  were  also  carried  out 
with  the  chamber  for  which  2f^  =  f , ^  .  They  confirmed  the  results  dis¬ 

cussed  here. 

Finally,  of  great  significance,  it  was  not  possible  to  achieve  self¬ 
sustenance  in  the  chamber  for  which  2f  ^  was  not  a  resonant  frequency. 

C_.  _  Theo re t i c a  1  Ana  lysis  of  the  Anemometer  Output  S i qua  1 

To  fully  understand  the  non  1 i near- feedback  experimental  results,  it 
was  necessary  co  analyze  the  anemometer  output  which  resulted  upon  exposure 
of  the  hot-wire  sensor  to  various  sound  fields. 

Considering  only  the  IT  mode,  the  bridge  voltage  of  the  anemometer 

was  calculated.  In  Figure  A  the  AC  bridge  voltage  is  plotted  versus 

dimensionless  time.  (The  curve  is  flat  in  the  vicinity  of  ,ot  =  0.  i  and 

2 n  becasue  of  the  assumption  that  whenever  the  calculated  values  of  E, 

using  Equations  (1).  are  less  than  E  ^ ,  the  no-flow  anemometer  output,  E, 

is  set  equal  to  E ^ .  This  means  that  during  periods  of  low  or  zero  air 

velocity,  the  bridge  voltage  is  equal  to  E  .)  Performing  a  Fourier  series 

o 

analysis,  the  harmonic  components  of  the  analytically  predicted  AC  bridge 
voltage  resulting  from  a  IT  sound  field  can  be  calculated.  They  are  shown 
in  Figure  5.  The  AC  signal  is  clearly  dominated  by  the  second  harmonic. 
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which  is  typic.il  of  the  rectification  process.  The  presence  of  the  first 

harmonic  is  primarily  due  to  (in*  temperature  difference  terms  [T  -  Tl  in 

w 

V 

Equat ions  (  1 ) .  * 

Considering  a  "distorted"  sound  field  produced  by  the  summation  of  a 
iT  and  a  -'T,  2l_  sound  field,  the  bridge  voltage  for  various  values  of  P- 
was  calculated,  where  P.<  i  is  the  ratio  of  the  amplitude  ot  the  sound  pres 
sure  of  the  T,  -L  sound  field  to  that  of  the  iT  sound  field  and  a  Fourie 
series  analysis  was  performed.  The  results  are  shown  in  Figure  6  and  7- 
It  can  be  seen  that  distortion  increases  the  component  of  E‘  at  the  funda 
mental  frequency  and  that  it  decreases  the  component  of  E1  .it  the  second 
ha  rmon i c . 


Figure  6.  Harmonic  Components  n  of  the  Analytically 

Predicted  AC  Bridge  Voltage  ol  t  lie  Anemometer 
for  P ■ i  =  0 .  0.5,  and  I  .  0 
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The  acoustics  of  the  chamber  are  important  here.  It  2t  is  a  resonant 
frequency  of  the  chamber,  a  relatively  strong  2f  acoustic  field  will  develop. 

If  it  is  not,  even  if  the  driving  at  this  frequency  is  strong,  the  acoustic 
field  at  2f  will  be  weak.  Also,  the  higher  harmonic  components  of  the  feed¬ 
back  signal  will  not  result  in  strong  acoustic  disturbances  at  these  fre¬ 
quencies,  first,  because  it  is  unlikely  that  they  are  also  resonant  fre¬ 
quencies,  and  second,  because  even  if  they  were  resonant,  it  is  relatively 
more  difficult  to  drive  the  higher  resonant  frequencies  than  the  lower 
resonant  f requeue i es . 

When  the  two  sound  fields,  the  initial  at  f  and  the  feedback  at  2f, 
combine,  the  resulting  sound  field  contains  both  frequencies.  The  hot-wire 
senses  this  new  sound  field  and  responds  with  an  AC  anemometer  voltage 
that  has  both  fundamental  and  second  harmonic  components.  Again,  the 
acoustics  of  the  chamber  play  an  important  role.  The  chambei  selectively 
responds  at  I  and  2f,  when  these  are  resonant  frequencies,  and  does  not  re¬ 
spond  well  to  the  higher  harmonic  components  of  the  input  signal.  The 
dependence  of  the  components  of  the  anemometer  signal  at  f  and  21  on  the 
relative  magnitude  of  the  acoustic  field  in  the  chamber  at  f  and  2f  was 
shown  in  Figure  7.  Figure  7  indicates  that  adding  even  a  relatively  lixv 
level  acoustic  field  at  2f  compared  to  the  initial  f  acoustic  field  results 
in  a  significant  component  of  the  feedback  signal  being  at  the  fundamental 
frequency  f. 

Again,  the  new  feedback  sound  field  is  superimposed  upon  the  existing 
sound  field  which  also  consists  ol  the  fundamental  and  second  harmonic  com¬ 
ponents,  and  a  new  feedback  signal  is  sent  to  the  driver.  Similar  existing* 
feedback  interaction  (c losed- loop)  processes  will  continue  as  long  as  the 
sound  field  exists. 

The  above  discussion  explains  why  an  acoustic  field  having  components 
at  both  f  and  2f  results  I rom  the  anemometer  analog  system.  The  two- frequency 
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character i st ic  is  a  result  of  the  nonlinear  nature  of  the  feedback  mechanism 


and  the  acoustics  of  the  chamber.  In  addition  to  the  factors  mentioned  above, 
the  amount  of  time  delay  or  phase  change  in  the  loop  is  also  of  importance. 

Tfie  more  out  of  phase  the  feedback  signal  is  with  respect  to  the  initial 
perturbation  in  the  chamber,  the  more  ampl i f icat ion  will  be  needed  for  self¬ 
sustenance  as  was  illustiated  by  Figure  3. 

Because  of  the  nonlinear  character  of  the  hot-wire  response,  all  the 
se 1 f -sus ta i ned  closed-loop  sound  fields  of  the  hot-wire  closed-loop  operation 
must  contain  higher  harmonic  components.  It  the  chamber  is  not  responsive 
at  the  2f  frequency,  i.e.,  if  2t  is  not  a  resonant  frequency  (or  close  to  a 
resonant  frequency) ,  the  sound  field  wi I  1  have  difficulty  being  sustained; 
hence,  the  geometry  of  the  chamber  is  or  primary  importance. 

More  extensive  discussions  of  the  experimental  closed-loop  investigations 
can  be  found  in  [8[,  f  1 2  f .  f  1  3 1  •  MM.  and  MM. 

Gasdynami ca 1 1 y  Induced  Second  Harmonic  Sound  Fields 

The  closed-loop  experimental  studies  indicated  that  the  second  harmonic 
sound  field  played  an  important  role  in  instability.  The  second  harmonic 
could  originate  from  two  souices.  It  could  be  directly  driven  by  the  pri¬ 
mary  energy  source  (combust  ion  in  rocket  engines;  the  acoustic  driver  of  the 
analog  process);  or  it  could  result  from  nonlinear  gasdynamic  effects  which 
i  arise  i  r'sonant  acoustic  tields  at  a  given  frequency  to  induce  sound  fields 
at  higher  harmonic  frequencies.  In  the  closed-loop  work  it  was  not  entirely 
cleat  how  much  second  harmonic  was  contributed  by  each  source;  hence,  an 
invest igat ion  was  made  of  gasdynami ca 1 1 y  induced  second  harmonic  sound  fields. 

The  object  of  these  experiments  was  to  determine  the  intensity  and  spatial 
distribution  of  t hr'  gasdynamical ly  induced  second  harmonic  sound  fields. 

Three  chamber  conf i gut  at  ions  were  used:  one  in  which  2f  _  -  f  ^  ^  ;  one 
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in  which  2f1^.  =  f  ^  ^  and  one  in  which  the  doubled  frequency  was  not  a 
resonant  frequency,  21^  =  f^.  Three  investigations  were  carried  out. 
First,  each  chamber  configuration  was  driven  at  f  ^  at  a  particular  sound 


f.'s),  at  five  different  sites,  to  determine  the  variation  in  amplitude  and 
phase.  Then,  one  location  was  chosen,  and  the  amplitude  of  the  iT  sound  field 
was  varied  and  the  level  of  the  result inq  2s  was  measured.  Lastly,  the 
location  of  the  driver  was  changed  to  determine  the  effect  of  driver  location. 

Initially  it  was  expected  that,  for  those  chambers  for  which  f2s  was  a 
resonant  frequency,  the  2s  would  have  the  amplitude  and  phase  characteristics 
of  the  resonant  sound  field.  This  was  found  to  be  only  partially  true.  The 
results  of  the  fiist  experiments  showed  the  gasdynamic  (2s)  sound  fields  were 
dependent  on  chamber  geometry-changing  the  length  of  the  chamber  changed  the 
spatial  distribution  of  the  second  harmonic.  Some  effort  was  directed  at 
characterising  the  various  2s  sound  fields.  The  distributions  appeared  to 
be  a  combination  of  the  resonant  sound  field  closest  in  frequency  to  f2S, 
plus  an  additional  component.  The  additional  component  appeared  to  be  in 
part  composed  of  a  radial  acoustic  mode  and  in  part  some  unknown  component. 
Also,  the  results  were  compared  to  the  analytically  predicted  gasdynamic 
sound  fields  obtained  by  Maslen  and  Moore  [  1 6 1 .  The  comparison  indicated 
the  2s  was  not  entirely  described  by  Maslen  and  Moore,  but  was  closer  to 
being  a  summation  of  the  Maslen  and  Moore  prediction,  plus  an  additional 
component  having  the  form  of  the  resonant  sound  field  closest  in  frequency  to 
f2s.  The  work  involved  in  trying  to  obtain  analytical  expressions  which 
would  characterise  the  measured  2s  sound  fields  was  extensive  when  compared 
to  the  usefulness  of  the  results  obtained,  so  these  efforts  were  terminated. 

The  results  of  the  second  experiments  showed  that  the  level  of  the  2s 
was  dependent  on  the  level  of  the  fundamental  and  that  to  obtain  a  significant 
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amount  of  nasdynamical ly  uenoratod  sci'oiul  hat  non  it  ,  an  intonse  fundamental 

sound  tit*  Id  would  In*  reuuired.  Typical  results  toi  the  t  t 

t  1  tl  ,  < l 

chamber  are  shown  in  Figure  8.  This  fioure  plots  the  ratio  ot  t  hi*  jntensi 
s'*  t  the  s  sound  field,  at  a  particular  lot.it  ion  in  (tie  shaT'ei,  to  the 
intensity  ot  t  ho  II  sound  field  onanist  t  tie  i  el  at  is  e  intensits  ot  the  II 
sound  t  i  e  I  si ,  when  compared  to  ambient  pressure.  There  were  on  I  s  relatisel 
small  variations  in  t  tie  results  foi  t  hi*  other  two  shamhei  t on t  i our  a t  i on s . 


I  ioure  8.  Relative  Amplitude  s't  (tie  iTasdvnamis  Sesond  Haimonit 
Sound  Field  as  a  Function  o  t  t  ho  Amplitude  s't  t  ho 
Fundamental  Sound  Field  at  Probe  Site  at  R  1.0  l  I  ’]. 
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When  the  location  of  the  driver  was  changed,  it  was  found  that  the 

overall  amplitude  of  the  second  harmonic  sound  fields  were  of  the  same 

\ 

order  of  magnitude  as  had  existed  previously.  A  variation  did  occur  in 
the  spatial  distribution  of  the  second  harmonic  in  some  cases.  The  cause 
of  this  is  unknown  but  is  probably  a  combination  ot  several  factors.  For 
instance,  there  were  always  slight  irregularities  in  the  driver  sound 
field--the  first  tangential.  These  irregularities  in  the  fundamental 
varied  with  driver  location  and  from  experiment  to  experiment.  Since  the 
fundamental  was  the  origin  of  the  gasdynamic  second  harmonic,  some  changes 
in  the  second  harmonic  were  to  be  expected.  Also,  the  driver  local  ion 
used  initially  was  a  good  location  for  actually  driving  a  second  harmonic 
resonant  sound  field;  the  later  driver  location  was  not. 

Further  details  of  the  second  harmonic  studies  can  be  found  in  [  1 7 1  . 

[ 1 8 1 ,  and  [19], 

The  above  studies  of  the  gasdynam i ca 1  I y  induced  second  harmonic  sound 
fields  gave  some  insight  into  the  nature  of  these  sound  fields,  but  did 
not  evaluate  the  contribution  of  these  sound  fields  to  instability.  A 
d  i  f f e ren t  .type  experiment  was  needed  for  that  purpose.  The  twin-chambers 
investigations  served  that  need. 

Tw in -Chambers  I nves t i gat  ions 

Simulation  of  actual  combustion  is  attained  in  the  expel imental  system 
by  operating  i  n  the  closed-loop  mode-- the  sensoi  responds  t  »■>  oscillations 
in  the  chamber  and  an  electrical  output  results,  which  is  amplified, 
phase-shifted  and/or  time  delayed,  and  the  resulting  signal  sent  to  an 
acoustic  driver  mounted  on  the  chamber,  which  in  turn  generates  new  osc i  I - 
lations,  and  so  on.  Altering  the  various  parameters  in  the  loop,  such  as 
amplification  or  phase,  alters  the  operating  charac t or i s t i c s ;  foi  example, 
the  oscillations  may  die  out  or  become  very  intense.  This  usually  occurs 


so  rapidly  that  it  becomes  difficult  to  make  detailed  studies  of  the 
phenomena  and  their  relationships  to  one  another.  To  overcome  these 
difficulties,  use  of  a  twin-chambers  system  was  initiated. 

Referring  to  Fiqures  9  and  10,  for  the  twin-chambers  studies,  the 
driver  on  the  first  chamber,  chamber  A,  was  activated  by  a  signal  origi¬ 
nating  from  some  combination  of  signals  from  two  function  generators. 

Ihis  resulted  in  oscillations  in  chamber  A,  to  which  the  sensor  responded. 
The  resulting  signal  was  then  amplified  and  sent  to  the  driver  on  the 
second  (or  twin)  chamber,  chamber  R.  which  resulted  in  oscillations  in  B. 
This  would  have  been  the  feedback  sound  field  in  A  if  chamber  A  were 
operating  closed-loop.  The  signals  to  t  lie  drivers  and/or  the  sound  fields 
in  each  chamber  were  then  compared  and  these  comparisons  indicated  whether 
or  not  a  single  chamber  operating  closed-loop  under  the  same  conditions 
would  be  stable  or  unstable.  (If  the  fundamental  and  second  harmonic 


Figure  9-  Schematic  Diagram  of  the  System  Used  foi  Pressure 
Sensitive  Open-Loop  Simulation  of  Closed-Loop 
Ope  rat  ion  [17]. 
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Figure  10.  Schematic  Diagram  of  the  Sys:em  Used  for  Velocity 
Sensitive  Open-Loop  Simulation  of  Closed-Loop 
Opera t i on  [  1 7 ] - 

frequency  components  of  the  feedback  sound  field  were  in  phase  with,  and 
had  the  same  or  greater  amplitudes  than  the  original  sound  field,  oscil¬ 
lations  would  be  sustained  or  magnified  if  the  system  were  operated 
c losed- loop. ) 

The  object  of  the  studies  was  to  examine  the  effect  of  second  harmonic 

content  on  the  operating  characteristics  of  the  two  types  of  systems--the 

one  having  a  linear  feedback  mechanism  (a  microphone  which  responds  linearly 

to  pressure  fluctuations)  and  the  one  having  a  nonlinear  feedback  mechanism 

(the  constant-temperature  hot-wire  anemometer  which  responds  nonlinearly  to 

velocity  fluctuations)--,  to  examine  the  importance  of  second  harmonic  to 

self-sustained  closed-loop  operation  and  to  determine  the  intensity  of  the 

second  harmonic  for  those  situations  for  which  self-sustenance  was  predicted. 

Three  sets  of  twin-chambers  were  used--the  2f  T  =  f  ,  the  2f,T  =  f,T 

1  I  *  I  ,  L  1  I  1  I  *  3  L 

and  the  2f  =  f  . 
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Fvii  tht-  I  inear  system,  referring  to  Figure  9,  the  procedure  was  ,is 

follows.  Initially  the  two  function  generators  were  set  to  produce  signals 

at  t  f  T  and  2t.  These  signals  were  summed  and  sent  to  the  driver  on 
1 1 

chamber  A.  The  microphone  sensed  the  resulting  sound  field  and  the  output 
of  the  microphone  was  amplified  and  sent  to  the  driver  on  chamber  B.  The 
phase  between  the  initial  signals  at  f  and  2f  were  adjusted  until  the  phase 
between  the  f  and  2t  components  ot  the  sound  field  in  A  was  the  same  as 
the  phase  between  f  and  2f  in  B.  Next,  the  amplifiers  were  adjusted  so  that 
the  amplitude  of  the  pressure  at  1  in  A  and  B  was  the  same.  The  amplitude 
of  the  signal  at  2f  to  the  summer  was  then  varied  between  ?ero  and  the  upper 
limit  of  operation  of  the  system  and  the  levels  of  the  resulting  sound 
fields  in  each  chamber  was  noted.  (After  this  adjustment  there  was  still  a 
relative  time  delay  between  the  sound  fields  in  the  two  v hamher s- - t he  time 
delay  necessary  in  the  loop  for  sustained  closed  loop  operation.) 

Figures  11  and  12  are  typical  of  the  results  obtained  for  the  linear 
system.  The  following  observations  can  be  made.  Referring  to  Figure  II, 
with  no  second  harmonic  voltage  to  the  driver  of  chamber  A,  primarily 
gasdynamic  second  harmonic  content  existed  in  each  chamber.  As  second 
harmonic  voltage  was  added  to  chamber  A,  a  lesser  amount  of  second  harmonic 
resulted  in  chamber  B.  The  amount  of  second  harmonic  in  chamber  B  was 
always  less  than  that  in  A.  Adding  second  harmonic  to  A  did  not  alter  the 
amplitude  of  be  fundamental  sound  field  i n  B.  The  only  conditions  under 
which  the  sound  fields  in  each  chambei  matched  occurred  when  A  was  not  being 
driven  at  the  second  harmonic  f reguency --when  primarily  gasdynamic  second 
harmonic  existed.  Figure  12  shows  the  amount  of  second  harmonic  content 
under  these  conditions  for  the  three  chamber  si.-es. 
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Figure  12.  Relationship  Between  Fundamental  and  Second  Harmonic 

Sound  Fields  as  a  Function  of  Intensity  of  Fundamental  [ 1 8 ] . 
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With  .1  lineai  system,  there  was  no  natural  mechanism  which  would 
enable  any  significant  second  harmonic  oscillation  othei  than  the  gas- 
dynamic  one  to  be  sustained  during  >.  losed-loop  operation.  It  a  per¬ 
turbation  occurred  in  the  chamber  at  the  second  harmonic  frequency  ,  this 
perturbation  would  tend  to  die  out.  The  second  hai  ’ionic  content  would 
not  be  necessary  tor  sel t -sustained  closed-loop  operation.  These  general 
results  did  not  strongly  depend  on  chamber  geometry;  similar  results  were 
obtained  tor  all  three  chamber  lengths.  The  results  are  consistent  with 
the  results  of  the  c  1  osed  -  I  oop  e\per  inrents. 

For  the  nonlinear  system,  referring  to  Figure  10,  the  procedure  was 
as  follows.  Initially,  the  two  function  generators  were  set  to  produce 
signals  at  f  =  f ,  and  2f.  These  signals  were  summed  and  sent  to  the 
driver  on  chamber  A.  The  hot-wire  sensed  the  sound  field  generated  and 
the  resulting  output  of  the  anemometer  was  amplified  and  sent  to  the  driver 
on  character  B.  The  phase  between  the  initial  signals  at  f  and  2f  was 
adjusted  until  the  phase  between  f  and  2f  of  the  sound  field  in  A  was  the 
same  as  the  phase  between  f  and  2f  in  B.  Next .  the  amplifiers  were  adjusted 
so  that  the  amplitudes  of  the  pressure  at  f  and  2f  in  A  matched  the  amplitudes 
of  the  pressure  at  f  and  2f  in  B,  respectively.  The  amplitude  of  the  signal 
at  2f  to  the  summer  was  then  varied  t  r'oin  ~ero  to  the  upper  limit  of  operation 
of  the  system  and  the  variation  in  the  results  sound  fields  recorded.  (Again, 
a  time-delay  existed  between  the  sound  fields.' 

Figures  12  and  13  are  typcial  of  the  results  obtained  tor  the  nonlinear 
system.  The  following  results  were  observed.  With  no  second  ha r. 'tonic  voltage 
to  the  driver  of  chamber  A,  the  sound  field  in  B  was  dominated  hv  second 
harmonic  content.  As  second  harmonic  voltage  was  added  to  A,  the  amount  of 
second  harmonic  content  in  the  sound  field  in  B  decreased  in  the  chambers 
for  which  2f,-|.  =  f ,  ^  at  2  t  ,  ^  =  f  y  'l  and  remained  approximately  the 
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Figure  13-  Levels  of  the  Sound  Fields  at  f  and  2f  in  Chambers  A  and  B 
as  a  Function  of  the  Voltage  at  2f  to  the  Acoustic  Driver 
on  Chamber  A;  Nonlinear  System;  2f,_  =  f  T  .  [ 1 7 ] * 

same  in  the  chamber  for  which  2f  =  fx.  There  was  an  amount  of  second 
harmonic  voltage  which,  when  added  to  A,  resulted  in  the  same  amplitude 
for  the  second  harmonic  sound  fields  in  both  chambers.  Adding  second 
harmonic  content  to  A  increases  the  amplitude  of  the  fundamental  sound 
field  in  B.  It  was  necessary  to  have  a  strong  second  harmonic  input  to 
chamber  A  in  order  to  obtain  a  strong  fundamental  sound  field  in  B.  There 
was  a  condition  under  which  the  sound  fields  in  each  chamber  had  the  same 
amplitudes  at  both  the  fundamental  and  second  harmonic  frequencies  and  the 
same  phases  between  fundamental  and  second  harmonic.  The  amount  of  second 
harmonic  existing  at  this  state  is  shown  in  Figure  12  for  each  of  the 
chambers.  The  results  for  the  2f;^.  =  tx  chamber  tended  to  be  somewhat 
errat ic. 
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With  t ho  nonlinear  system,  a  sound  field  at  t ho  fundamental  t ret uency 
naturally  resulted  in  a  response  primarily  at  the  harmonic  frequency.  The 
system  would  not  bo  able  to  sustain  itself  under  these  condit  ions.  It 
would  bo  necessary  to  have  a  driver  input  at  both  the  fundamental  and 
second  harmonic  frequencies  in  order  for  the  sound  fields  in  chambers  A 
and  B  to  match.  Second  harmonic  content  would  be  necessary  for  so  I f-susta i nod 
closed-loop  operation.  This  was  consistent  with  the  results  of  the  closed- 
loop  investigations. 

What  was  not  entirely  consistent  with  the  closed-loop  results  was  the 
prediction  that  the  chamber  for  which  2f  was  not  a  resonant  frequency 
could  be  self-sustaining.  With  the  twin-chambers  experiments  it  is  predicted 
that  a  2  f ;  .j.  =  fx  chamber  could  be  self-sustaining  il  enough  feedback  ampli¬ 
fication  were  available  in  the  loop.  (A  factor  which  is  not  noted  in  the 
figures  is  that  it  took  approximately  twice  as  much  power  to  achieve  the 
"matched"  situation  for  the  2f^  =  fx  chamber  as  it  did  for  either  the  ?f,j  = 

f  _  or  2f  _  =  f  _  ,  chambers.)  The  twin-chambers  results  indicate 

l  T  ,  T  t  T  T  ,  ?  L 

that,  if  the  closed-loop  experiments  had  had  a  more  powerful  amplifier  and 
an  acoustic  driver  able  to  withstand  higher  power  input,  perhaps  closed- 
loop,  se I f-susta i ned  operation  would  have  been  achievable. 

Further  discussion  of  the  twin-chambers  investigations  can  be  found 
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This  ri'si’.ii  i  h  was  begun  because  of  .1  desire  to  dptetmine  whothor  01 
not  o'ft.i  in  effects  th.it  hud  boon  observed  expel  imentally  could  bo  pro* 
dieted  mat hom.it  i c .1 1  1  y .  It  was  toll  that  this  would  ho  helpful  in  ex¬ 
plaining  the  expet  imental  results  to  other  engineers  and  scientists  and 
in  i nterpretat ion  ot  the  observed  phenomena.  In  ordot  that  the  pie- 
dictions  of  the  analytical  work,  could  he  easily  understood,  it  was  desired 
to  avoid  methods  which  requited  t he  s i mu  1 1 aneous  numei  ical  solution  of  a 
system  ot  partial  differential  equations.  lot  this  to.tson,  the  method 
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It  was  decided  to  use  a  similai  methodology  foi  the  analysis  ot  velocity- 
sens  it  ive  ciwthus tic'ti  instability  problems. 

It  was  tound  that  if  a  set  ot  approx  ini.it  ions  slightly  different 
from  those  used  in  [ 1 0 1  was  employed,  a  wave  equation  could  he  dot  ived 
which  contained  gasdynamic  nonlineai ities  ot  all  oidets,  rather  than  only 
t hose  of  order  two  as  in  [ 1 0 ] .  The  assumptions  involved  in  this  analysis 
are  discussed  by  Peddieson,  Ventt  ice,  and  Pin  dy  1  1  i' I  and  tuither  details 
are  given  by  Wong  [?0l.  The  appropriate  wave  equation  is  (in  d imens i on  1 es' 
t  orni) 

•*„  ’  *  «  '  '  »  1  «  •  •'  -  "  -  "• 

(.1  J>  ♦  •  C 4“ / 2 >  1  t  W  [I  -  (y  -  ♦  \\-  •  \\x  ,M|1'  '  ^  ''  0 

where  V  is  the  gradient  operator,  T  is  the  velocity  potential,  W  is  the 
burning  rate,  y  is  the  ratio  ot  specific  heats  foi  a  pet  tec t  gas  and  t 


is  t  i  me  . 


To  be  consistent  with  t  tie  wot  k  ot  previous  investigators,  and  to 


simplify  the  computation  eltort,  it  w.is  decided  to  develop  .1  second-older 
form  of  (4).  Ttiis  was  done  by  assuming  that 

T  -  1  IM')  *■  ■f  ( x .  y  ,  .  Ol 


(5) 


W  =  w  ( 7 I  r  i  w  ( x  ,  y  ,  r  ,  t  1 


where  1  is  a  measure  ('I  the  initial  d  i  st  urbnnces ,  7  is  the  .1xi.1l  coordinate 
(parallel  to  the  main  flow)  in  the  combustor,  and  the  guant  it  ies  with  the 
superposed  bar  are  associated  with  (lie  mean  How.  Substituting  (5)  into 
(4),  as sumirig  w  =  0(e)  and  w  -  0(i  ),  neglecting  all  terms  ot  0(t  ’)  and 
dividing  by  1  ,  one  obtains 

T  a  +  w.yj>  1  2uT  4-  -  V>  i  ,  (2Vf  •  C .1 1  +■ 


where 


(1  -  1)V‘\*>.T  +  w)  -  0 


u  (z)  ■  /  wdl 


(6) 


(7) 


is  the  steady-state  gas  velocity.  The  pertui bat  ions  of  velocity,  density 
pressure,  and  temperature  associated  with  ^  are  respectively 

u  =  V  -I- 


(8) 


p  -  -  ( ,1(  «J>  1  U,1  J  +  t  (ti  -  (7  -  >  )  (  vf1)  ■  )  /?  1 

p  -  >  l  T  ^  *  u  .'t4*  ■*  r  (u  -  (;'t«f)‘)/2l 
T  -  -  ( ^  -  11  ( T  j  T  +  uT^  r  tu  /.’) 

It  was  first  attempted  tv'  use  t  lit'  vapot  i  rat  ion  limited  combustion  model 
of  Pt  iem  [l]  in  conjunction  with  (6).  In  the  present  notation  this  can  be 
wr  i  1 1  en  as 

w  -  ( w  / 1  ^  1(1  -  T  ( f  / T )  (^l  1  -  T  ) '  1 


(h) 
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where  u^  is  the  velocity  of  the  liquid  drops.  When  (9)  is  combined  with 
(6),  the  result  is  a  single  equation  for  if .  Because  of' the  non-polynomial 
non  I i nea r i t  i  es  in  (9),  this  equation  cannot  be  analyzed  conveniently  by 
the  Galerkin  method.  For  this  reason,  another  version  of  the  method  of 
weighted  residuals,  the  orthogonal  collocation  method,  was  used.  The 
derails  of  this  method  are  discussed  in  [18]  and  [20].  It  was  eventually 
concluded  that  this  method  was  unreliable  when  applied  to  this  type  of 
problem  and  was  abandoned.  It  is  believed  that  this  unreliability  is 
related  to  the  fact  that  the  nonlinear  terms  in  (6)  have  the  correct  form 
only  for  moderate  values  of  t: .  This  can  be  illustrated  by  the  following 
example. 

Consider  a  nonlinear  elastic  system  governed  by  the  equation 

f  +  f  +  ,  N  (f.  f)  =  0  (10) 

where  f  is  the  displacement,  N  is  a  nonlinear  function,  i  <<  1 ,  and  a 

superposed  dot  indicates  differentiation  with  respect  to  time.  Suppose 

N  ( f  ,  f)  =  N  ff  +  0(,  )  (II) 

o 

Then  the  second  order  form  of  (10)  would  be 

f  +  f  +  tN  ff  =  0  (12) 

o 

This  equation  can  be  integrated  once  to  yield 

log [ ( 1  +  i N  f)/(l  +  cN  f  )  ]  -  i  N  (f  -  f  )  -  r2  N  2  (f2  -  f  2)/2  =  0  (13) 

o  o  o  o  o  o  o 

where 

f (o)  =  f  ,  f (o)  =  f  (O) 

o  o 

For  cN^  <<  1,  (13)  can  be  simplified  to 

f2  +  f2  -  2.  N  f 3 /3  +  0  ( r  2  N  2)  =  f  2  +  f  2  (IS) 

O  OOO 


This  is  consistent  with  the  original  approximat i on  leading  to  (12).  For 

iN  <<  1  it  is  expected  that  the  solution  of  (12)  will  be  a  resonable 

\ 

approximat ion  to  the  solution  of* (10).  Equation  (15)  indicates  that  the 
solutions  are  periodic  as  one  would  expect  for  free  vibrations  of  an 
elastic  system. 

Suppose  that  is  sufficiently  large  so  that  the  inequality  cN^  <<  1 
is  not  satisfied.  Then  (12)  is  not  a  rational  approximation  to  (10)  and 
the  solution  of  (12)  may  not  have  the  characteristics  of  an  elastic  system. 
As  an  example,  consider  cN^  >>  1.  Then  (13)  can  be  simplified  to 

f  =  f  exp  [e2N2(f2-f2)/2]  (16) 

o  o  o 

For  f  =0  this  becomes 
o 

f  =  f  exp (e2  N  2  f 2 /2)  (17) 

o  o 

In  this  case  the  solution  is  not  periodic  and  grows  without  bound.  The 
application  of  the  collocation  method  to  (6)  produces  a  set  of  nonlinear 
ordinary  differential  equations  governing  the  values  of  <j>  at  the  collocation 
points.  These  equations  contain  nonlinear  terms  similar  to  the  last  term 
in  (12).  It  appears  that  most  choices  of  collocation  points  cause  the  co¬ 
efficients  of  some  of  these  terms  (equivalent  to  >:N  )  to  become  sufficiently 

o 

large  so  that  the  equations  predict  behavior  similar  to  that  indicated  by 
(17)-  This  happens  even  in  the  absence  of  combustion  terms.  It  may  be 
that  more  success  could  be  achieved  by  applying  the  collocation  method 
directly  to  (4),  but  this  matter  was  not  pursued. 

Sin,.1  the  collocation  method  could  not  be  made  to  produce  satisfactory 
results,  it  was  decided  to  employ  the  Galerkin  method.  As  mentioned  pre¬ 
viously,  this  method  cannot  be  used  in  conjunction  with  (9).  Thus,  (9) 
was  replaced  for  purely  velocity  sensitive  combustion  with 


w 


w  n  u 


(18) 


where  n  is  colled  the  interaction  index.  This  is  analogous  to  the 
customary  treatment  of  pressure  sensitivity  (see,  for  instance,  Powell  [10]). 

As  in  the  case  of  pressure  sensitivity,  (18;  can  be  generalized  to  account 
for  history  effects  as 

w*wn(u‘_u‘)  (19) 

where  u  (x,  y,  z,  t)  =  u(x,  y,  z,  t  -  t)  and  x  is  called  the  time  lag. 

Equation  ( 1 8)  can  be  related  to  burning-rate  laws  meant  to  apply  to  special 
types  of  combustion  processes  such  as  (9)-  for  an  example  of  this,  see 
Wong  [20,  21].  With  the  use  of  (19).  it  was  possible  to  achieve  significant 
progress. 

Peddieson,  Wong,  and  Ventrice  [19]  discussed  the  application  of  (6) 
and  (19)  to  transverse  wave  motion  in  a  cylindrical  combustion  chamber. 

An  approximate  solution,  based  on  the  IT,  1R,  and  2T  acoustic  modes  of 
the  chamber,  was  used  and  the  Galerkin  method  was  employed  to  carry  out 
a  mod a  1  analysis.  This  resulted  in  a  system  of  ordinary  differential  equa¬ 
tions  that  could  be  solved  to  determine  the  temporal  behavior  of  the  modal 
amplitudes.  This  was  done  numerically  using  a  fourth-order  Runga-Kutta 
method.  (For  further  details,  see  [ 1 9l  and  [20,  21].)  Stability  boundaries 
were  determined  and  the  behavior  of  the  pressure  at  the  chamber  wall  was 
determined.  Some  typical  stability  boundaries  are  shown  in  Figures  1 4 - 1 7 . 

It  can  be  seen  that  all  boundaries  resemble  rectangular  hyperbolas  in  the 
n-e  plane.  A  purely  velocity-sensitive  system  will  always  be  stable  to 
infinitesimal  disturbances.  Thus,  the  instabilities  associated  with  Figures 
14-17  ore  triggered  instabilities. 

It  was  felt  that  the  numerical  solutions  discussed  above  could  be  usefully 
supplemented  by  using  the  method  of  multiple  scales  (see, 
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for  instance, 
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Nayfeh  [22))  to  obtain  approximate  analytical  solutions  to  the  governing 
equations  for  the  modal  amplitudes.  To  i  nves t  i  ga  te -t  he  feasibility  of 
doing  this,  Googerdy  [23.  2^]  considered  purely  pressure-sensitive  com¬ 
bustion  and  replaced  (19)  by 

w“-wn(3^-3  <frt)/e  (20) 

Rather  than  deal  with  a  cylindrical  chamber,  he  considered  the  mat hema t i ca 1 1 y 
simpler  problems  of  transverse  motion  in  a  narrow  rectangular  chamber  and 

I 

transverse  motion  in  an  annular  chamber  with  a  narrow  gap  width.  The  first 
problem  is  analogous  to  that  of  radial  oscillations  in  a  cylinder,  while 
the  second  is  analogous  to  the  problem  of  combined  radial  and  transverse 
motion  in  a  cylinder.  Both  of  these  cases  were  treated  numerically  by 
Powell  [10].  Using  the  method  of  multiple  scales  it  was  possible  to  produce 
most  of  Powell's  results  in  closed  form  and  to  provide  physical  explanations 
for  them.  This  gave  confidence  to  the  use  of  the  perturbation  method. 

McDonald  [25,  26]  then  applied  the  method  of  multiple  scales  to  the 
problem  of  purely  velocity-sensitive  combustion  instability  in  an  annulus 
of  narrow  gap  width.  Transverse  motions  were  treated.  For  standing  waves, 
the  simplest  solution  able  to  show  the  effect  of  quadratic  non  1 i near i t i es 
has  the  form 

v)>  ((' ,  t )  =  f  i  ( t )  cos  0  +  f  2  ( t )  cos  20  (2  0 

where  0  is  a  polar  coordinate  which  loactes  radial  cross  sections  in  the 

annulus.  Substituting  (21)  and  r=l  into  (6)  and  (19).  using  the  Galerkin 

method  to  carry  out  the  modal  analysis,  and  assuming  that  the  gasdynamic 

nonlinearities  are  negligible  compared  to  the  combustion  non  1 i near i t ies , 

results  in  the  equations 

f  \  +  f\  +  e  i  ( cr  f  i  +  2wn  f  i  f  ’ )  *  0 

(22) 

f;  +•  f ,  +  c(afj  -  wn  f  x '  /2 )  =  0 


^7 


where  =  w/  .  Solving  these  equations  by  the  method  of  multiple  scales, 
subject  to  the  initial  conditions 

f,(0)  =  f.-(0)  =  6,  ti(  0)  =  1,  f  -  (0>  =  0  (23) 

leads  to  the  un i form! y-va I i d  first  approx ima t i ons 

f'l  =  exp(  -  w  t/2)  sec  fen  (1  -  exp(  -  w  t/2))/2  '  ]  sin  t  +  0(e) 

(2*0 

f;  =  -  (1/2  /  )  exp  (  -  w  t/2)  tan  [rn(l  -  exp(  -  wt/2))/2  ^  ]  ■ 
sin  2 1  +  0(c) 

These  simple  solutions  illustrate  the  basic  features  of  the  system.  Both 
the  secant  and  the  tangent  become  infinite  when  their  arguments  take  on 
the  value  -t/2.  Thus,  if  en  <  tt  2  '*  ,  the  arguments  of  the  secant  and  tangent 
in  (2*0  are  less  than  -t/2  for  all  t,  and  fi  and  t;  approach  zero  as  t  ap¬ 
proaches  infinity.  I  f  in  %  rr  2  ,  on  the  other  hand,  these  arguments  reach 

t/2  in  a  finite  time,  causing  fi  and  f;  to  become  infinite  in  a  finite  time. 
The  former  situation  is  stable,  while  the  latter  is  unstable.  The  equation 
of  the  stability  boundary  in  the  n-c  plane  is 

l  /: 

n  =  2  t/c  (25) 

For  this  situation  it  can  be  shown  that 

lim  fj  =  (2/t)  sin  t 

t  v  CO 

1  (26) 
lim  f :  =  -  ( 1  /2  ^  !t)  sin  2t 

£  -*■  co 

Thus,  a  neutral  oscillation  is  approached  when  mass  addition  due  to  unsteady 
combustion  is  exactly  balanced  by  mass  loss  due  to  steady-state  flow. 
McDonald  [25,  26]  has  found  similar  solutions  for  traveling  waves  and  show 
that  the  conclusions  discussed  above  are  qualitatively  unchanged  by  the 
inclusion  of  gasdynamic  nonlinearities.  He  has  further  shown  that  the 
qualitative  aspects  of  the  perturbation  solutions  are  in  agreement  with 
those  of  direct  numerical  calculations  based  on  (22).  This  work  illustrates 


vividly  the  usefulness  of  the  combustion  law  ( 1 8 ) .  Closed-form  solutions 
similar  to  (24)  could  never  be  obtained  using  a  combustion  law  of  the 
form  (9) . 
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Cone ! us i ons 


The  research  has  resulted  in  useful  experimental  analog  and  analytical 
techniques  for  the  study  of  Reynolds  number  dependent  processes.  The 
techniques  are  general  and  could  be  applied  to  any  Reynolds  number  dependent 
process.  During  the  course  of  the  research,  application  to  combustion  in¬ 
stability  in  liquid  propellant  rocket  engines  was  emphasized. 

Experimentally,  the  initial  object  of  the  work  was  to  study  the  use¬ 
fulness  of  a  constant- temperature  hot-wire  anemometer  as  an  analog  tool 
for  the  i nves t i ga t ion  of  combust  ion  instability.  To  accomplish  this, 
analytical  and  experimental  open-loop  studies  of  the  anemometer  were  carried 
out  and  compared  to  analytical  open-loop  studies  by  Heidmann  (2,  3]  of 
vaporization  limited  combust  ion. (Open- loop  studies  compare  existing  pressure 
perturbations  and  the  resulting  perturbations  in  the  burning  rate  in  the 
case  of  combustion,  or  in  the  anemometer  output  in  the  case  of  the  analog, 
in  order  to  predict  actual  closed-loop  operation.)  The  comparisons  indicated 
qualitative  agreement  among  the  three  studies--the  second  harmonic  frequency 
distort ial  of  the  fundamental  frequency  perturbation  was  the  important 
factor  in  the  development  of  instability. 

In  closed-loop  studies,  the  burning  rate  in  combustion,  or  anemometer 
'utput  in  the  analog  process,  feeds  back  upon  itself,  resulting  in  new 
pressure  perturbations,  which  combine  with  the  previous  perturbat ions 
to  form  a  new  acoustic  field,  and  so  on.  The  analog  was  refined  to  run 
closed-loop  using  either  the  original  nonlinear  analog  device--the  anemometer 
or  a  linear  analog  device--a  microphone.  Experiments  were  then  run  to 
investigate  the  validity  of  the  open-loop  predictions  of  closed-loop  behavior 
These  studies  confirmed  the  importance  of  second  harmonic  content  to  insta¬ 
bility  in  the  nonlinear  case  and  gave  insight  into  the  reasons  for  this 
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importance-- second  harmonic  content  is  a  natural  result  of  the  nonlinearity 
of  the  process  and  is  necessary  to  sustaining  perturbations  at  both  the 

X 

fundamental  and  second  harmonic  frequencies. 

Two  important  similarities  in  the  charac ter  i  s t ics  of  combustion  in¬ 
stability  and  the  closed-loop  analog  process  were  identified.  first,  in 
combustion  instability,  a  minimum  perturbation  amplitude  is  necessary  to 
excite  instability.  In  the  analog  this  was  also  true.  Second,  in  com¬ 
bustion  instability,  an  initial  sinusoidal  perturbation  becomes  harmonically 
distorted  once  instability  develops.  Again,  in  the  analog  this  was  also 
found  to  be  true.  These  studies  also  indicated  that  chamber  resonance  might 
play  a  more  significant  role  in  causing  instability  than  had  been  anticipated. 
Second  harmonic  content  was  not  found  to  be  important  in  the  linear  feedback 
process . 

The  closed-loop  analog  studies  did  not  fully  indicate  the  difference 
in  character  between  linear  (pressure  sensitive)  and  nonlinear  (velocity 
or  Reynolds  number  sensitive)  processes,  so  a  twin-chambers  technique  was 
developed  for  this  purpose.  In  the  twin-chambers  work,  instead  of  having 
the  analog  device  feeding  back  on  itself,  the  feedback  signal  was  sent  to 
an  identical  twin  chamber  so  that  the  originating  and  feedback  perturbations 
could  be  studied  separately.  These  studies  indicated,  even  more  clearly 
than  the  closed-loop  studies,  that  second  harmonic  was  an  integral  part  of 
a  nonlinear  process  and  necessary  to  the  sustenance  of  the  fundamental 
frequency  perturbations.  They  also  indicated  that  this  was  not  true  for 
the  linear  case--the  linear  case  was  primarily  a  single  frequency  process. 

All  of  the  experimental  results  indicated  that  the  nonlinear  aspects 
of  the  combustion  process  were  the  dominant  mechanisms  control inq  combustion 
i  n  s  t  a  b  i  I  i  t  y . 
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Both  the  closed- loop  and  twin-chambers  analog  technique  of  studying 
instability  proved  to  be  valuable  methods  for  gaining  Thsight  into  the 
phenomena  being  studied  and  both  have  potential  for  application  to  other 
s i m i 1 ar  problems . 

The  analytical  work  established  a  method  for  computing  stability 
boundaries  for  velocity-sensitive  combustion.  It  was  established  that  the 
method  used  moderate  amounts  of  computer  time  and  predicted  results  that 
were  in  qualitative  agreement  with  the  known  characteristics  of  velocity- 
sensitive  combustion  instability.  The  results  of  the  numerical  work  were 
substantiated  by  comparison  of  approximate  solutions  obtained  by  perturbation 
procedures.  In  some  cases,  the  perturbation  solutions  were  obtained  in 
closed  form.  It  is  felt  that  these  explicit  solutions  are  very  helpful 
in  understanding  the  physical  nature  of  velocity-sensitive  combustion 
instabi 1 i ty. 
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All  references  marked  with  *  are  papers  or  reports  which 
resulted  from  this  research.  In  addition  to  the  reports  listed, 
semiannual  reports  were  submitted  every  six  months.  Two  of  these 
semiannual  reports,  July  through  December  1973  and  January  through 
July  1975,  were  of  special  significance  and  were,  therefore,  as¬ 
signed  CR  numbers.  They  are  listed  in  the  references  as  numbers 
6  and  12.  Also,  informal  reports  were  submitted  on  a  somewhat 
monthly  basis  to  the  project's  Technical  Monitor,  Dr.  Richard  J. 
Priem.  These  reports  have  some  details  in  them  that  do  not  appear 
in  the  published  work. 
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